Resonance-assisted decay of nondispersive wave packets 
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We present a quantitative semiclassical theory for the decay of nondispersive electronic wave 
packets in driven, ionizing Rydberg systems. Statistically robust quantities are extracted combin- 
ing resonance assisted tunneling with subsequent transport across chaotic phase space and a final 
ionization step. 
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The accurate and quantitative modeling of open quan- 
tum systems at high spectral densities remains a serious 
challenge for quantum theory and computational physics. 
Even the ab initio treatment of apparently simple, atomic 
one or two electron systems under nonperturbative forc- 
ing, in the energy range of laboratory experiments, has 
become available only during the last two decades 0. 
Furthermore, it appears a safe bet that the driven three 
body Coulomb problem Q will be the most complicated 
"many particle" quantum problem still amenable to a 
(numerically) exact and complete solution, even on the 
most powerful supercomputers. Thus, alternative theo- 
retical strategies which allow a precise treatment, with- 
out saturating our computational resources nor loosing 
quantitative predictive power, are in need. 

Specifically in the range of high spectral densities, 
semiclassical approaches - which try to deduce the spec- 
tral structure of a given quantum system from the un- 
derlying Hamiltonian dynamics - open such an avenue. 
However, while semiclassics of bounded systems is well 
developed , its generalization for open systems Q has 
not been fully accomplished yet: Our semiclassical under- 
standing of tunneling and decay phenomena in systems 
with mixed regular-chaotic classical phase space struc- 
ture remains rather incomplete, and semiclassics still has 
to prove its potential to come up with robust (a priori) 
quantitative predictions for specific experiments. 

In the present Letter, we improve on that situation, by 
elaborating a fully quantitative semiclassical treatment of 
a paradigmatic example of open system dynamics at high 
spectral densities, in the context of light-matter interac- 
tion and coherent control: Nondispersive wave packets 
have been shown to be ubiquitous in periodically driven 
quantum systems with underlying mixed regular-chaotic 
phase space 0, IE 13 ■ They can be launched along 
essentially arbitrary phase space trajectories, which can 
be manipulated in real time 9\. In contrast to "tradi- 
tional" nondispersive wave packets, which rely on the 
harmonicity of the spectra they are built on |lfj, these 
strongly localized eigenstates of the driven system are 
robust against perturbations 0, 0, EI> what is a con- 
sequence of the Kolmogorov-Arnold-Moser theorem 
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FIG. 1: Classical phase space of driven hydrogen, gener- 
ated by the Hamiltonian Q, for u>o = 1 and Fo — 0.041, 
and spanned by the action-angle variables of the unperturbed 
atom, at cot — 0. Nondispersive quantum wave packets are 
localized at the center of the principal elliptic island (con- 
centric structure in the middle of the plot) and move along 
the resonantly driven Kepler orbit sketched in the inset (the 
arrow shows the direction of the driving force at the outer 
turning point), in phase with the driving field. For the cho- 
sen parameters, a prominent, higher order nonlinear 5 : 1 
resonance strongly affects the tunneling coupling of the prin- 
cipal island's eigenstates to the surrounding chaotic sea. This 
resonance is expressed by the five-fold structure right within 
the outermost torus which confines the principal island. 



The characteristic property which defines their unique 
potential for quantum control purposes (including quan- 
tum memory applications [l2|. i.e. .within the framework 
of molecular quantum computing 13]) is their essentially 
eternal life time on experimentally relevant time scales 
0, El- The life time is determined by the decay rate 
of the wave packet from the elliptic island to which it is 
anchored to in classical phase space (see Fig.^l, i.e., it is 
essentially limited by a finite quantum mechanical tun- 
neling rate thr oug h classically impenetrable phase space 
barriers 0,QIIB|- 

However, while always very long, the precise value of 
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the wave packet's life time is not a statistically robust 
quantity - it fluctuates over several orders of magnitude 
under tiny changes of the boundary conditions, which 
may be stray fields or other uncontrolled imperfections in 
the experiment. This sensitivity stems from the coupling 
of the wave packet to the chaotic component of phase 
space, into which it is completely embedded (see Fig.^l. 
We will derive a quantitative scmiclassical estimate for 
the mean value and the variance of the wave packet's 
decay rate, without adjustable parameters. 

For the sake of clarity, we consider a simplified (yet 
reliable 7] ) formulation of the problem we have in mind: 
Atomic hydrogen initially prepared in a highly eccen- 
tric (extremal parabolic) Rydberg state aligned along the 
field polarization axis of a linearly polarized microwave 
driving field of amplitude F and frequency uj can be de- 
scribed by the Hamiltonian 



H(p z ,z) 
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+ Fz cos(wi) , 



(1) 



where the configuration space of the electron has been 
confined to the z-axis (for the actual three dimensional 
problem, such confinement can be achieved by applying 
a weak additional static electric field along the polariza- 
tion axis of the periodic drive 0). The classical phase 
space structure generated by this Hamiltonian is illus- 
trated in Fig. ^ an d is completely fixed by the choice of 
the scaled driving field frequency loq — ujuq and ampli- 
tude Fq — Ftlq, respectively, as a consequence of the scale 
invariance of Hamilton's equations of motion. Specifi- 
cally, we set wo = 1 and Fq = 0.041, in order to realize 
the essential ingredients of our subsequent treatment: (i) 
an elliptic island (at the center of Fig. [I}, which is (ii) 
completely embedded into the chaotic sea, and (iii) en- 
closes a prominent, higher order r : s nonlinear resonance 
right within its outermost confining torus. 

For different values of no, and accordingly adjusted u> 
and F, we find (by exact numerical diagonalization of 
the Floquet Hamiltonian defined by 0) JJ) wave packet 
eigenstates localized on the elliptic island, together with 
their decay rates T wp . Under changes of no, the latter 
(which, in this specific physical realization, are nothing 
but ionization rates mediated by multiphoton coupling 
to the atomic continuum) fluctuate wildly, as shown in 
Fig. H below. 

To determine the mean value and the variance of F wp , 
let us start from the random matrix ansatz for an effective 
Hamiltonian employed in (ljj : 
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The matrix structure comprises the wave packet state 
localized at the center of the regular island (with en- 



ergy E wp ), and N additional states located in the chaotic 
phase space region, one of which being subject to decay 
with the rate r c . The chaotic states are assumed to be 
strongly coupled to each other, such that the correspond- 
ing sub-block (Hij ) can be described by a random matrix 
from the Gaussian orthogonal ensemble [l^. In accor- 
dance with the resonance-assisted coupling mechanism 
to be explained below, we assume that the wave packet 
is dominantly coupled to only one of those chaotic states, 
via the (comparatively small) matrix element V fc . As a 
qualitative, conceptional improvement over |l5| . where 
V TC and r c were free fitting parameters, we here derive 
quantitative estimates for these two coupling constants. 

Within our semiclassical framework, transport to the 
continuum is only possible by traversing a large part 
of chaotic phase space. Quantum mechanically, this is 
equivalent to multiphoton transition amplitudes which 
couple the initial state - the wave packet - to Rydberg 
states with ionization potentials smaller than u>, i.e., with 
principal quantum numbers n above n c = n 3 / 2 / \/2loq. 
The final, ionizing step occurs from these highly excited 
states by a one-photon process. In the effective Hamilto- 
nian of (0), this final ionization process is accounted for 
by attributing a finite decay rate r c to one of the chaotic 
states, which is well estimated by the Golden Rule ex- 
pression 



T r = 0.265F 2 w" 10/3 n7 3 ~ 1.26 x 10~ 3 r 
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a.u. 



(3) 



for the single-photon ionization rate of a Rydberg state 
with quantum number n c [l7| . 

A random-matrix average over the eigenvalues and 
eigenvectors of the chaos block (Hij ) gives rise to the 
Cauchy-type probability density 
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of the wave packet's decay rate |l5|. The characteristic 
scale To is derived as To = (irV rc /uj) 2 r c , exploiting that 
the chaotic eigenvalues are uniformly distributed within 
an energy interval of width to, due to the inherent pe- 
riodicity of the Floquet spectrum. Because of the large 
fluctuations implied by the distribution the relevant 
statistical quantity is given by the averaged logarithm 
|l8| . An explicit calculation of this geometric mean gives 



exp[(ln(r wp ))] = r = 



r c 



(5) 



which depends only on the single-photon ionization rate 
T c from ©, on the driving frequency u, and on the cou- 
pling matrix element V IC between the regular island and 
its chaotic surrounding. V rc can be derived using purely 
semiclassical arguments, as follows. 

Let us recall that nonlinear resonances between the 
external drive and the local modes of a regular phase 



3 



space region induce higher order perturbative couplings 
between the locally quantized eigenstates of that region 
[l9|. This leads to the phenomenon of resonance- assisted 
tunneling which, originally proposed for near-integrable 
dynamics, can be generalized to mixed regular-chaotic 
systems as well |2Cj. There it provides the dominant 
semiclassical mechanism for the tunneling process that 
connects the "ground state" of a regular island (localized 
at the island's center) to states within the chaotic sea 
(which, in turn, are coupled to the atomic continuum 
with a rate r c , given in Eq. ©). 

Quantitatively, the classical dynamics near a r:s reso- 
nance (s oscillations match r driving periods) is described 
by the effective pendulum Hamiltonian |l9l | 



H 1C 
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+ 2V r:s cos(rd) , 



(6) 



expressed in terms of the local action-angle variables /, 9 
of the principal elliptic island, where s is absorbed in the 
definition of 9. The central action I r - S of the r:s reso- 
nance, and the effective mass and coupling parameters 
w r:s , V r:s , can be extracted with little numerical effort 
from the classical phase space, by computing the area 
covered by the separatrices of the r:s resonance, and by 
evaluating the stability matrix of the associated periodic 
points j20]. Quantizing ©, we see that the perturba- 
tion term induced by the r:s resonance couples the A:th 
excited state of the principal island, given by the plane 
wave (9\k) ~ exp(ik9) in angle space, to the states \k±r), 
with a strength V r - S . Since the wave packet state corre- 
sponds to the principal island's ground state with k = 0, 
its coupling to the chaotic sea is thus given by the effec- 
tive matrix clement 

V^vJflj^, (7) 

1=1 ^0 — &br 



with E k = (7 fe - I r:s ) 2 /(2m r:s ) and T k = h(k + 1/2). 
Here, \l C 'f) denotes the first state within this perturbative 
sequence of higher order couplings that is located outside 
the principal island, and the prefactor V r:s accounts for 
the final step from |(Z C — 1) • r) to \l c ■ r). 

In our exemplary case of Fig. ^ the driving induces 
a prominent 5:1 resonance within the principal regular 
island on which the wave packet is localized. From the 
area covered by the principal island - which is determined 
numerically - we deduce no ~ 135 as the critical quan- 
tum number at which the island's fifth excited state (i.e., 
the first state to which the ground state is coupled by the 
5:1 resonance) is located exactly on the outermost invari- 
ant torus. For no < 135, we can therefore identify this 
fifth excited state with the first basis state of the chaos 
block in the Hamiltonian J5J), and set V rc = V r - S , which, 
together with Eqs. (|3I5|I and the numerically computed 
value of V r:s , leads to our central result 
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FIG. 2: (color online). Comparison of the exact quantum de- 
cay rates of the wave packet states (symbols) with the semi- 
classical prediction for their mean value r wp [solid line, given 
by Eqs. JHJ and J5J] and the standard deviation [dashed lines, 
see Eq. JT0J]. The rates are plotted as a function of the princi- 
pal quantum number no of the Rydberg orbit along which the 
wave packet is launched. The quantum results were obtained 
for the same field parameters as Fig. Open circles and di- 
amonds distinguish different basis sizes form the numerical 
diagonalization. For no < 150, the basis was chosen such as 
to resolve all bound states with an ionization potential larger 
and equal to the one-photon energy w. For no > 150, only 
bound states below 2no were numerically resolved. 



for the mean decay rate of the wave packet. At larger 
principal quantum numbers (no > 135), the coupling 
from the regular island to the chaotic sea involves one 
more perturbative step in Eq. Q, which substantially 
reduces the matrix element V rc - This leads to 



T wp ~ 1.7 x 10 
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(n - 131) 



r a.u. , no > 135 . 



(9) 



As a consequence, the semiclassically calculated decay 
rates exhibit a sudden drop at no ~ 135. 10 describes 
the monotonic decrease (mind the expression's denomi- 
nator!) of the rates with no up to no ~ 250, where an 
additional drop is expected to occur by virtue of yet an- 
other perturbative step in Eq. (JJJ. 

Fig. [3 highlights that our semiclassical result T wp fits 
the exact quantum decay rates T wp of the wave packet 
very well. On average, after a first exponential decrease 
for no < 60, the decay rate settles to a plateau at 
T wp ~ 10~ 13 a.u., and eventually drops again around 
no > 135. Beyond that, our model also succeeds to very 
well characterize the quantum rates' fluctuations by their 
logarithmic standard deviation derived from fllf: 



r. 



9.6 x 10 



" 13 n 1 ^ 2 a.u. , no < 135, 



^/varpn^wp)] - J ( [ln(r wp ) - ln(r wp )] 2 ) = n . (10) 

This universal result is indicated by the dashed lines in 
(8) Fig. |21 which are given by r wp e ±7r and which encompass 



4 



the overwhelming part of the data. Despite the sensitive 
dependence of individual decay rates on boundary condi- 
tions and classical phase space structures, our combined 
semiclassical and random matrix approach provides a ro- 
bust statistical characterization, over a wide range of ex- 
citation energies. 

The disagreement between our semiclassical result and 
the quantum data below uq ~ 60 is not surprising, since, 
at such low excitations, the principal island cannot ac- 
commodate a sufficient number of quantum eigenstates to 
resolve the 5 : 1 resonance (approx. r/2 states are needed 
for satisfactory resolution |20j). Thus, a necessary con- 
dition (item (iii) above) for the resonance assisted tun- 
neling mechanism to manifest itself is not satisfied at too 
small no < 60. 

In conclusion, the plateau-like structure born out in 
Fig. 21 in the range n ~ 60 . . . 135, corroborates the rel- 
evance of the resonance assisted tunneling mechanism for 
the decay properties of nondispersive wave packets, in an 
experimentally routinely accessible energy range 0, 0|. 

— 1/2 

Note that the average trend T wp ~ n , predicted by 
(JSJ and confirmed in Fig. |3 is markedly different from 
an on average exponential decay with uq, as assumed in 
[l5ll2l| . Indeed, a similar step structure arises in earlier 
published data on nondispersive hydrogen wave packets 
|21| (at a slightly different field amplitude), which, on 
the basis of our presently improved understanding, is also 
attributed to resonance-assisted tunneling. Thus, while 
the precise extension of the plateau depends on the spe- 
cific choice of parameter values, and certainly also on 
the effective dimensionality of the electron's configura- 
tion space, the structure as such is a robust fingerprint 
of resonance assisted tunneling, as well as the universal 
value 7r of the decay rates' standard deviation. An ex- 
perimental verification of this prediction is challenging, 
though appears in reach for the most advanced experi- 
mental setups to date @, 0] • 
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